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Abstract 

The AI = 1/2 rule and direct CP violation e'/e in kaon decays are studied within the framework 
of the effective Hamiltonian approach in conjunction with generalized factorization for hadronic ma- 
trix elements. We identify two principal sources responsible for the enhancement of Re^o/R'e^2: 
the vertex-type as well as penguin-type corrections to the matrix elements of four-quark operators, 
which render the physical amplitude renormalization-scale and -scheme independent, and the non- 
factorized effect due to soft-gluon exchange, which is needed to suppress the AI = 3/2 K ^ tttt 
amplitude. Contrary to the chiral approach which is limited to light meson decays and fails to 
reproduce the A2 amplitude, the aforementioned approach for dealing with scheme and scale issues 
is applicable to heavy meson decays. We obtain ReAo/ReA2 = 13 — 15 if ms(lGeV) lies in the 
range (125—175) MeV. The bag parameters Bi, which are often employed to parametrize the scale 
and scheme dependence of hadronic matrix elements, are calculated in two different renormalization 
scehemes. It is found that and Bq^\ both of order 1.5 at = 1 GeV, are nearly 75 scheme 
independent, whereas B^}^ 7 as well as B^f^ show a sizable scheme dependence. Moreover, only 
Bi'l 4 exhibit a significant dependence, while the other i?-parameters are almost m<j indepen- 
dent. For direct CP violation, we obtain e'/e = (0.7 - 1.1) x lO'^ if ms(lGeV) = 150 MeV and 
e'/e = (1.0 - 1.6) X IQ-^ if 

nis is as small as indicated by some recent lattice calculations. 
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I. INTRODUCTION 



The celebrated A J = 1/2 rule in kaon decays still remains an enigma after the first 
observation more than four decades ago. The tantalizing puzzle is the problem of how 
to enhance the A0/A2 ratio of the AI = 1/2 to AI = 3/2 K vrvr amplitudes from 
the outrageously small value 0.9 [see Eq. ( [5.1D below] to the observed value 22.2 ± 0.1 
(for a review of the AI = 1/2 rule, see [||]). Within the framework of the effective weak 
Hamiltonian in conjunction with the factorization approach for hadronic matrix elements, 
the Aq/A2 ratio is at most of order 8 even after the nonfactorized soft-gluon effects are 
included |jl| . Moreover, the fi dependence of hadronic matrix elements is not addressed in the 
conventional calculation. In the past ten years or so, most efforts are devoted to computing 
the matrix elements to 0{p^) in chiral expansion. This scenario has the advantages that 
chiral loops introduces a scale dependence for hadronic matrix elements and that meson loop 
contributions to the Aq amplitude are large enough to accommodate the data. However, 
this approach also exists a fundamental problem, namely the long-distance evolution of 
meson loop contributions can only be extended to the scale of order 600 MeV, whereas the 
perturbative evaluation of Wilson coefficients cannot be reliably evolved down to the scale 
below 1 GeV. The conventional practice of matching chiral loop corrections to hadronic 
matrix elements with Wilson coefficient functions at the scale n = (0.6 — 1.0) GeV requires 
chiral perturbation theory and/or perturbative QCD be pushed into the regions beyond their 
applicability. 

Another serious difficulty with the chiral approach is that although the inclusion of chiral 
loops will make a large enhancement for Aq, it cannot explain the A2 amplitude. For example, 
in the analysis of in which a physical cutoff Ac is introduced to regularize the quadratic and 
logarithmic divergence of the long-distance chiral loop corrections to K ^ vrvr amplitudes, 
the amplitude A2 is predicted to be highly unstable relative to the cutoff scale Ac and it even 
changes sign at Ac > 650 MeV . In the approach in which the dimensional regularization 
is applied to regularize the chiral loop divergences and to consistently match the logarithmic 
scale dependence of Wilson coefficients, the predicted A2 amplitude is too large compared 
to experiment ^ , indicating the necessity of incorporating nonfactorized effects to suppress 
the A J = 3/2 amplitude This implies that not all the long-distance nonfactorized 
contributions to hadronic matrix elements are fully accounted for by chiral loops. In short, 
it is not possible to reproduce Aq and A2 amplitudes simultaneously by chiral loops alone. 

Even if the scale dependence of K —>■ matrix elements can be furnished by meson 
loops, it is clear that this approach based on chiral perturbation theory is not applicable to 
heavy meson decays. Therefore, it is strongly desirable to describe the nonleptonic decays 
of kaons and heavy mesons within the same framework. 

In the effective Hamiltonian approach, the renormalization scale and scheme dependence 
of Wilson coefficients is compensated by that of the hadronic matrix elements of four-quark 
operators 0{n) renormalized at the scale /x. Since there is no first-principles evaluation 
of (0(/i)) except for lattice calculations, it becomes necessary to compute the vertex- and 
penguin-type corrections to (O) (not {0{fi)) !), which account for the scale and scheme de- 
pendence of {0{fi)), and then apply other methods to calculate (O). The AI = 1/2 rule 
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arises from the cumulative effects of tlie sliort- distance Wilson coefficients, penguin opera- 
tors, final-state interactions, nonfactorized effects due to soft-gluon exchange, and radiative 
corrections to the matrix elements of four-quark operators. As shown in Q, the last two 
effects are the main ingredients for the large enhancement of Aq with respect to A2. 

Contrary to the nonfactorized effects in charmless B decays, which are dominated by 
hard gluon exchange in the heavy quark limit and expected to be small due to the large 
energy released in the decay process, the nonfactorized term in K tttt is anticipated to 
be large and nonperturbative in nature, namely it comes mainly from soft gluon exchange. 
One can use tt+tt^ to extract the nonfactorizable contributions to the hadronic matrix 

elements of (V — A){V — A) four-quark operators @]. 

Instead of using scheme- and scale-independent effective Wilson coefficients, one can al- 
ternatively parametrize the hadronic matrix elements in terms of the bag parameters B^^^ and 
B\ which describe the scale and scheme dependence of hadronic matrix elements {Qi{fj.)). 
These non-perturbative parameters are evaluated in the present paper. We have checked 
explicitly that these two seemingly different approaches yield the same results. 

The prediction of e' /e in the standard model is often plagued by the difficulty that the 
result depends on the choice of the renormalization scheme. Presumably this is not an issue 
in the effective Wilson coefficient approach. Unfortunately, as we shall see in Sec. IV, our 
predictions for e'/e are scheme dependent for reasons not clear to us. 

The present paper is organized as follows. In Sec. II we construct scheme and scale 
independent effective Wilson coefficients relevant to K ^ ttti decays and direct CP violation 
e'/e. The bag parameters Bi are evaluated in Sec. III. Based on the effective Wilson 
coefficients or bag parameters, K titt amplitudes and direct CP violation are calculated 
in Sec. IV and their results are discussed in Sec. V. Sec. VI is for the conclusion. 



II. EFFECTIVE WILSON COEFFICIENTS 

The effective Hamiltonian relevant to K ^ tttt transition is 

Kfr(A5 = 1) = Ec.(/^)Q.(/^)1 +h-c., (2.1) 

where 

= + r?/i(/x), (2.2) 

with r = -VuV;JiVudV:,), and 

Q3(5) = isd)y_^J2(9Q)v-A{v+A)y Qm = isadp)y_j,^{qf3qa)v-Aiv+A)^ (2-3) 

9 Q 

3 3 

^7(9) = 9 (^^)v-A (^qi.m)v+AW-A)^ Qsm = 7^{sadp)y_^ ^ eg(g^g„) ^^^(^_^) , 

with Q^-Qq being the QCD penguin operators, Qj-Qw the electroweak penguin operators 
and {(liq2)v±A = ^i7m(-'- ^ 1^)12- The sum in Eq. ( p.3| ) is over light flavors, q = u,d, s. 
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In the absence of first-principles calculations for hadronic matrix elements, it is custom- 
ary to evaluate the matrix elements under the factorization hypothesis so that {Q{fi)) is 
factorized into the product of two matrix elements of single currents, governed by decay 
constants and form factors. However, the naive factorized amplitude is not renormalization 
scale- and 75 scheme- independent as the scale and scheme dependences of Wilson coefficients 
are not compensated by that of the factorized hadronic matrix elements. In principle, the 
scale and scheme problems with naive factorization will not occur in the full amplitude since 
{Q{fi)) involves vertex-type and penguin-type corrections to the hadronic matrix elements 
of the 4-quark operator renormalized at the scale /x. Schematically, 

weak decay amplitude = naive factorization + vertex— type corrections (2-4) 

-|- penguin— type corrections + spectator contributions + ■ ■ ■ , 

where the spectator contributions take into account the gluonic interactions between the 
spectator quark of the kaon and the outgoing light meson. The perturbative part of vertex- 
type and penguin-type corrections will render the decay amplitude scale and scheme in- 
dependent. Generally speaking, the Wilson coefficient c(/i) takes into account the physics 
evolved from the scale Mw down to /i, while {Q{^)) involves evolution from fj, down to the 
infrared scale. Formally, one can write 

{Q{li))=g{lx,fif){QM), (2.5) 

where /ij is a factorization scale, and g{fi,fif) is an evolution factor running from the scale 
fi to fif which is calculable because the infrared structure of the amplitude is absorbed into 
{Qil^f)). Writing 

c^'^if^f) =c{fi)g{fi,fif), (2.6) 

the effective Wilson coefficients will be scheme and /x-scale independent. Of course, it appears 
that the /i-scale problem with naive factorization is traded in by the /i/-scale problem. 
Nevertheless, once the factorization scale at which we apply the factorization approximation 
to matrix elements is fixed, the physical amplitude is independent of the choice of /x. More 
importantly, the effective Wilson coefficients are 75-scheme independent. In principle, one 
can work with any quark configuration, on-shell or off-shell, to compute the full amplitude. 
Note that if external quarks are off-shell and if the off-shell quark momentum is chosen 
as the infrared cutoff, g{fi,fif) will depend on the gauge of the gluon field p|. But this 
is not a problem at all since the gauge dependence belongs to the infrared structure of 
the wave function. However, if factorization is applied to {Q{fif)), the information of the 
gauge dependence characterized by the wave function will be lost. Hence, as stressed in 
PJTO[], in order to apply factorization to matrix elements and in the meantime avoid the 



gauge problem connected with effective Wilson coefficients, one must work in the on-shell 
scheme to obtain gauge invariant and infrared finite c^^ and then applies factorization to 
(<5(/i/)) afterwards. Of course, physics should he fif independent. In the formalism of 
the perturbative QCD factorization theorem, the nonperturbative meson wave functions are 
specified with the dependence on the factorization scale fif 0. These wave functions are 
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universal for all decay processes involving the same mesons. Hence, a consistent evaluation of 
hadronic matrix elements will eventually resort to the above-mentioned meson wave functions 
determined at the scale /i/. 

In general, the scheme- and jU-scale-independent effective Wilson coefficients have the 



where /x j is the factorization scale arising from the dimensional regularization of infrared di- 
vergence , and the anomalous dimension matrix as well as the constant matrix fy arise 
from the vertex-type corrections to four-quark operators. For kaon decays under considera- 
tion, there is no any heavy quark mass scale between mc and mx- Hence, the logarithmic 
term emerged in the vertex corrections to 4-quark operators is of the form ln(/ij//i) as shown 
in Eq. (|2.7|) . We will set yUj = 1 GeV in order to have a reliable estimate of perturbative 
effects on effective Wilson coefficients. 

It is known that the penguin operators Qs e do not induce tttt directly, but their 

Fierz transformations via {V — A){y ^ A) —2{S + P){S — P) do make contributions. 
Applying equations of motion, (Q^fiifJ,)) are proportional to mj^/[ms{fi) +mg(/i)] with q = u 
or d. This means that, contrary to current x current operators, the matrix elements (Q^^ei^)) 
for K — 7CTT transition under the vacuum insertion approximation do exhibit a fi dependence 
governed by light quark masses. The fi dependence of the Wilson coefficients c^fi{fi) is 
essentially compensated by that of light quark masses (the cancellation becomes exact in 
the large- A'^c limit). Of course, the near cancellation of /i dependence does not imply that 
factorization works for the matrix elements of density x density operators since the scheme 
dependence of c^fi^fi) still does not get compensation. It is thus advantageous to apply 
the aforementioned effective Wilson coefficients to avoid the scheme problem caused by 
factorization. And in the meantime, the fif dependence of C5^(/i/) is largely canceled by 
that of quark masses entering the matrix elements {Q5fi{fif))- 

To proceed, we note that the renormalization-scale and -scheme independent effective 
Wilson coefficient functions zf^ are given by (for details, see P,p!0[]): 



form 11,121: 



CM 




Cj{fi) + penguin— type corrections. 



(2.7) 
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If = zrifi) + 



;eff 



zeS 
^10 



Air 



as 
Air 



a. 



8i 



in 



a 



9i 



Sir 



{21 



Likewise, the effective Wilson coefficients have similar expressions in terms of except 
that ?/i = 2/2 = 0. The reason why we call the l.h.s. of Eq. ( p.8|) as z^^ rather than zf 
will become shortly. In Eq. (|2.8| ) the superscript T denotes a transpose of the matrix, the 
anomalous dimension matrix 7^ as well as the constant matrix fy arise from the vertex 
corrections to the operators Qi — Qio, Ct, Cp and Ce from the QCD penguin-type diagrams 



of the operators Qi,2, the QCD penguin-type diagrams of the operators 
electroweak penguin- type diagram of Qi,2, respectively: 
-2 



Qq, and the 



Ct 



a, 



-K 

3 
4 

-K 

3 



8 /2 

- -K 

9 V3 



G{mu,k,ij,)j zi, 
G{ms, k, fi) - G{md, k, jj)^ Z3 

G{mu,k,fi)] (2:1 + 3^2), 



E 

i=u,d,! 



G{mi){zi + ze) 



(2.9) 



where k is a parameter characterizing the 75-scheme dependence in dimensional regulariza- 
tion, i.e., 

1 NDR, 
HV, 

in the naive dimensional regularization (NDR) and 't Hooft-Veltman (HV) schemes for 75, 
and the function G{m, k, jj) is given by 



(2.10) 



dxxil — x) In 



^ — k'^xil — xY 



(2.11) 



with fc^ being the momentum squared carried by the virtual gluon. The explicit expression 
for 7y is given in |Tl|. For reader's convenience, we list here the constant matrix fy [p!0|,|6[|: 
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(2.12) 
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in the NDR scheme, and 
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(2.13) 



in the HV scheme. Note that the 66 and 
have been corrected in and ||13|| . 



entries of fy given in [|T^] are erroneous and 



The results of a direct calculation of z^^ and y^^ in NDR and HV schemes using Eq. 
are displayed in Table I. Formally, the effective Wilson coefficients are scale and scheme 



independent up to the order a^. This implies that the scheme independence of requires 
that the Wilson coefficients Zi{fi) appearing in the vertex- type corrections and in Ct, Cp and 
Ce be replaced by the lowest-order (LO) ones zf'~', while Zi appearing in Ct be the next- 
to-leading order (NLO) one, and likewise for y^^. Therefore, we define the effective Wilson 
coefficients z'^^ and y'^^ similar to Eq. ( |2.8| ) except for the above-mentioned replacement, for 
example. 



„eff 



Air 



(2.14) 



li 



From Table I we see that the scale independence of z'^^, which is good to the accuracy of the 
third digit, is significantly better than for example, the values of 2;4,6(NDR) and 2;4,6(HV) 
to NLO are quite different, but their effective Wilson coefficients zf^Q are obviously scheme 



independent. 



Note that zf, 



^eff 



are enhanced relative to their NLO values by about 



three times as they receive large corrections proportional to OsCt. By contrast, the scheme 
independence of |/^^ is not as good as z^^. In particular, yf^, which plays an important role 
in e'/e, shows a slight scheme dependence for reasons not clear to us. 



The effective Wilson coefficients appear in the factorizable decay amplitudes in the com- 



binations 



-2i 



+ l^^lf-i and a2i_i 



, J_^eff 
1 AT, 



5). For K ^ Tin decays 



nonfactorizable effects in hadronic matrix elements can be absorbed into the parameters a! 



eff 
i 



14-16 



^eff 
hi 



+ [j^ + X2i ] Z^^-l, 



"2i-l 



with Xi being the nonfactorizable terms. Likewise 

1 



"2i 



y2f + + X2^) y^U, 



(2.15) 



(2.16) 
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TABLE I. AS =1 Wilson coefficients at /x = 1 GeV for mt = 170 GeV and = 325 MeV, 

^ MS ' 

taken from Table XVIII of [17]. Also shown are the effective Wilson coefficients zf^ and zf^ (see 
the text), yf^ and y'-^ in NDR and HV schemes with /x = 1 GeV, /i/ = 1 GeV and fc^ = rn\/2. 
Note that yi = y2 = 0. 





TO 




XI V 


Zi 1,-LNiJllJ 


Zi V-n- V j 


Zi -LJK'j 


Zi l-ti V J 




1 zL'^'^ 

1 .^00 




1.0/1 


1 fil 4 


1 fi7S 
i.u 1 


1 71 8 
1 . 1 10 


1 71 
1 . / 10 




n 7/18 

-U. / 4:0 




-U.U4U 


-1 .uzy 


1 089 


1 1 1 
-1.110 


1 1 1 n 

- 1 . 1 lU 


^3 






007 


03Q 


034 


039 


039 


za 


-0.012 


-0.035 


-0.017 


-0.080 


-0.085 


-0.081 


-0.084 




0.004 


0.008 


0.004 


0.024 


0.024 


0.024 


0.025 




-0.013 


-0.035 


-0.014 


-0.094 


-0.086 


-0.086 


-0.086 


Zf/a 


0.008 


0.011 


-0.002 


0.025 


0.047 


0.063 


0.069 


zs/a 


0.001 


0.014 


0.010 


0.025 


0.016 


0.016 


0.013 


Zg/a 


0.008 


0.018 


0.005 


0.039 


0.057 


0.072 


0.078 


zio/a 


-0.001 


-0.008 


-0.010 


-0.015 


-0.012 


-0.011 


-0.012 




LO 


NDR 


HV 


yf (NDR) 


yf (HV) 


(NDR) 


yf (HV) 


2/3 


U.Uoo 


u.uoz 


U.Uo 1 




u.uou 


u.uou 


u.u^y 


2/4 


-0.061 


-0.058 


-0.061 


-0.051 


-0.055 


-0.053 


-0.053 


2/5 


0.013 


-0.001 


0.016 


0.004 


0.003 


0.003 


0.002 


2/6 


-0.113 


-0.111 


-0.097 


-0.161 


-0.130 


-0.160 


-0.138 




0.036 


-0.032 


-0.030 


-0.051 


-0.019 


-0.052 


-0.028 


2/8/a 


0.158 


0.173 


0.188 


0.287 


0.295 


0.285 


0.300 


2/9/a 


-1.585 


-1.576 


-1.577 


-2.013 


-1.919 


-2.053 


-1.948 


2/io/a 


0.800 


0.690 


0.699 


1.339 


1.205 


1.355 


1.216 



To proceed, we shall assume that nonfactorizable effects in the matrix elements of {V — 
A)(y + A) operators differ from that of {V — A){y — A) operators; that is, 

Xll = Xi = X2 = Xz = Xi = X<s = Xio, 

Xlr = Xh = X% = X'J = Xs, (2.17) 

and Xlr 7^ Xll- Theoretically, a primary reason is that the Fierz transformation of the 
(y — A){V + A) operators 05,6,7,8 is quite different from that of {V — A)(y — A) operators 
Oi, 2,3,4 and Og^io [|10|- Experimentally, we have learned from nonleptonic charmless B decays 
that Xlr{B) 7^ Xll{B) ||T0| , p!3| . As shown in 0, the nonfactorized term xll can be extracted 
from — s> TT+TT^ decay to be 

XLL = -0.73. (2.18) 

Contrary to the nonfactorized effects in hadronic charmless B decays, which are dominated 
by hard gluon exchange in the heavy quark limit ^ and expected to be small due to the 
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large energy released in the decay process, the nonfactorized term in — vrvr is large and 
nonperturbative in nature, namely it comes mainly from soft gluon exchange. 

III. NON-PERTURBATIVE PARAMETERS 5, 



In the literature it is often to parametrize the hadronic matrix elements in terms of 
the non-perturbative bag parameters and B^"^^ which describe the scale and scheme 
dependence of the hadronic matrix elements {Qi{fi)): 



(<5i(/^))o 



{Qi)o^^ {Qi)2 



VIA ' 



(3.1) 



where (Qi)"^^^ denote the matrix elements evaluated under the vacuum insertion approxima- 
tion. In order to evaluate the parameters Bf'''^\ as an example we consider the vertex- type 
and penguin-type corrections to the hadronic matrix element of the four-quark operator Qi 



in the NDR scheme fTO]: 



1 + 



a. 



in 



_2 In ^ + 3 



(Qi(/i/)) + ^f61n^-9j (Q2(^/)) 



(3.2) 



G(m„, k, fi) - (P(/i/)) - ^ (G{mu, /c, /x) - ^ ) (Qrifif) + QgifJ^f)), 



where 



P = ~Q3 + Q4-^Q5 + Q6- 



(3.3) 



The parameters B^^''^^ are then determined from Eq. (|3.2| ). In general, the hadronic param- 
eters Bf'''^\n, Hf) have the expressions:* 



(0) 



(Qi)o + ^ U In ^ + f-vj {Q^)o - g (G(mJ - ) (P)o 



li 



2 
-/ 
3 



(2) 



^(0,2) 



{Ql)2 + ^[lvl^^+fvj {Q^)2 

(G(m„) - {Q, + Q9)2|/(Qi)y^^, 
(Q2)o,2 + ^(7yln^ + f^) (Q.)0,2 



*Note that our convention for Qi,Q2 (and hence Bi,B2) differs from that in |T^,17| where the 
labels 1 and 2 are interchanged. 
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(G(m.) - (Q7 + Q9)o,2|/(Q2)S^, 



47r \ /i 



3i 



5f = |(Q4)o + ?^f7vln^ + ryl (g.)o 



47r y /i 



4j 



[G(m„) + G(m,) + GK)] (P)o /(g4)r, 



pf = I (Q5)o + £ (^v In ^ + ) (g.)o [/(g5)r 



5j 



Air \ /i 



6i 



= {{Qs)o,2 + ^(lv\n^ + fv] (g.)o4/(Q8)S^ 



8i 



= l{Qio)o,2 + ^(lv\n^ + fy] (g.)o.2|/(gio)S^. (3.4) 



47r y jj, 



Wi 



For simplicity we have dropped the parameters k and /i in the argument of the function G. 
Note that the effective Wilson coefficient z^^ is in general not equal to Zi{fi)Bi{^), but the 
physical amphtude in terms of z*^^ or Zi{fi)Bi{fi) is the same. 

The — > vrvr matrix elements under the vacuum insertion approximation read (see e.g., 

0) 



NrJ ' ' 3 \ N, 



' c 



3 V NJ' '^^'^ 3 \ 

(Q7)o = ^/i.^;^ + Ix, {Qr)2 = - ^X, 
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(Q8)o = VQJkv' + ^X, 



(Q8)2 = V3Ut 



X, 



(^9)0 

(Qio)o 



Ix (i - ± 



Nr 



(Oio)2 = -]=X (1 + ^ 



V2 



Nr 



(3.5) 



where X 



mz 



and 



2 



m„ + md md + rUs 



2 2 



(3.6) 



characterizes the quark-order parameter (qq) which breaks chiral symmetry spontaneously. 

To evaluate Bf'''^^ we need to take into account nonfactorized effects on hadronic matrix 
elements. As discussed in Sec. II, this amounts to replacing 1/Xc by 1/Xc + Xll for (V — 
A){y — A) quark operators and by 1/Nc + Xlr for {V — A){V + A) operators. For example, 



{Qi)o 



{Q2)o 



l-6x 



LL, 



(Q5)0 



l + 3x 



LR- 



(3.7) 



° to be 



Although the nonfactorized term xll is fixed by the measurement of vr ' vr 

—0.73 0, no constraint on xlr can be extracted from —>■ vrvr. Nevertheless, we learned 
from hadronic charmless B decays that xlr 7^ Xll As shown in Fig. 1, the parameters 



^5, and eJ^^ are quite sensitive to Xlr, whereas Bq, Bj'^ and Bj^^ stay stable. Lattice 



calculations suggest that B5 Bq = 1.0 ± 0.2 and B^^^ = 0.6 ± 0.1 at /i = 2 GeV P[. The 



lattice results roughly imply the constraint —0.45 < xlr < 0. We shall see in Sec. IV that 
the Aq/A2 ratio and e'/e are not very sensitive to the variation of Xlr- Using rriu = 3.5 
MeV, rud = 7.0 MeV, = 140 MeV at = 1 GeV and xlr = —0.1, the numerical values 
of R-^'^'* are listed in Tables II and III. 




Xlr 



FIG. 1. The bag parameters (solid line), Bj (dotted line) and Bj (dashed line) evaluated 
in the NDR scheme versus xlr for 771^(1 GeV) = 140 MeV. 
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TABLE II. Numerical values of the non-perturbative bag parameters at fj, = 1 GeV in 
NDR and HV schemes for m,(l GeV) = 140 MeV, /// = 1 GeV and = m\/2, xll = -0.73 and 
Xlr = —0.1. The results for sf'^ in the absence of nonfactorized contributions (i.e. xll = Xlr = 0) 
are shown in parentheses. 







sf 


r(0) 


4°' 


^(0) 




^(0) 


^9 


^10 


NDR 2.5(1.8) 


8.7(3.1) 


-0.1(2.^ 


]) 2.7(2.7) 


1.0(1.3) 


1.5(1.5) 


1.0(1.3) 


1.6(1.6) 


3.1(1.5) 


3.1(1.5) 


HV 2.6(2.1) 


8.0(2.7) 


0.9(3.8 


) 2.6(2.7) 


1.7(2.0) 


1.5(1.5) 


1.7(1.9) 


1.6(1.6) 


2.9(1.4) 


2.9(1.4) 



TABLE III. Same as Table II except for B\ 







^(2) 




o(2) 


^(2) 


^(2) 


NDR 


0.34(0.75) 


0.34(0.74) 


1.0(1.4) 


1.6(1.6) 


0.35(0.76) 


0.35(0.76) 


HV 


0.37(0.81) 


0.36(0.79) 


1.8(2.1) 


1.6(1.6) 


0.37(0.81) 


0.37(0.81) 



IV. K ^TTTT ISOSPIN AMPLITUDES AND e'/e 



In terms of the effective Wilson coefficients defined in Sec. II, the CP-even AI = 1/2 
and AI = 3/2 K ^ tttt amphtudes have the form 0: 



ReAo = 



Gf Re(KdK:. 



\/2 cos 5o 



2 1 1 

-ai - -a2 + 04 + -{aj - ag + aio) X 



ReAo = 



v/2 cos ^2 l-VL 



IB 



ai + 02 + 2^~(^7 + 09 + aio, 



^X + VSfyasl, (4.1) 



where = ^2^/^42 describes the isospin breaking contribution to K'^ n^n^ due to the 
n — 1] — 1]' mixing, and 6q as well as S2 are S-wave tttt scattering isospin phase shifts. For 
simplicity, we have dropped the superscript "eff" of aj in Eq. (|4.1| ). 

The direct CP- violation parameter e'/e given by the general expression 



e 

e 



can be recast in the form 



e^ _ Gfuj 
e ~ 2|e|ReAo 



10 /ImA2 ImAo 
^\e\ VReAa Re^o 



(64 + - + ^&10)X 



(4.2) 



[ cos (5o 

1 1 



ijJ COS ^2 



1 - ^^ib) 
67 + 69 + 610 )X/V2 + V3i;V^6s 



(4.3) 
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where uj = A-ij = 1/22.2 . 

Alternatively, the K ^ tctc amphtudes and direct CP violation can be expressed in terms 
of the non-perturbative parameters Bf''^^: 

V2 cosdo ~{ 

V2 cos 02 1 - i Im ~i 

and 

e 2|£:|Ke/io [cosoo^g 

--^j:yM'\Q^)T''\- (4.5) 

CO COS (32 J 

We have checked expHcitly that the numerical values of AI = 1/2, 3/2 amplitudes and 
e'/e obtained using the parameters sf'''^^ given in Tables II and III are the same as that 
calculated using the effective Wilson coefficients zf^ and y^^ shown in Table I, as it should 
be. Since z"^^ and yf^ show a better scale independence than z^^ and y^^, we will use the 
former set of effective Wilson coefficients in ensuing calculations.^ 

Using fiiB = 0.25 ±0.02 0, 5o = (34.2 ±2.2)°, ^2 = -(6.9 ±0.2)° |2|, xll = -0.73 and 
Xlr = —0.1, we plot in Fig. 2 the ratio A0/A2 as a function of at the renormalization 
scale n = I GeV. Specifically, we obtain 



ReAo 
ReA, 




1 GeV) = 125 MeV, 

1 GeV) = 150 MeV, (4.6) 
1 GeV) = 175 MeV. 



It is clear that the strange quark mass is favored to be smaller and that the prediction is 
renormalization scheme independent, as it should be. In Fig. 3 we study the dependence of 
Aq/A2 on the unknown nonfactorized term Xlr- It turns out that the ratio decreases slowly 
with Xlr, but it stays stable within the allowed region —0.45 < xlr < 0. 
For direct CP violation, we find for lm{VtdVtl) = 1.29 x 10"^ (see Fig. 4) 



e 
e 



1.56 (1.02) X 10-3 at (1 GeV) = 125 MeV, 

1.07 (0.70) X 10-=^ at m, (1 GeV) = 150 MeV, (4.7) 
0.78 (0.51) X 10-3 at (1 GeV) = 175 MeV, 



''^The results of Aq/A2 and e'/e calculated using the bag parameters B^'"^^ and effective Wilson 
coefficients zf^ and yf^ are numerically very similar except for e'/e in the HV scheme which is 
slightly small in terms of i?-parameters by around 15% compared to that evaluated in terms of yf^ . 
For example, Eq. (|]|) leads to e'/e = 0.59 x 10"^ at ms{l GeV) = 150 MeV in the HV scheme, 
while Eq. (|]l]) yields e'/e = 0.70 x lO'^ [see Eq. (|^] in the same scheme. 
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FIG. 2. The ratio of ReAo/Reyl2 versus (in units of GeV) at the renormahzation scale /U = 1 
GeV, where the sohd (dotted) curve is calculated in the NDR (HV) scheme and use of xlr = —0.1 
has been made. The solid thick line is the experimental value for ReAo/Re^2- 



< 15 " 
<U 

E 

O 10- 



-0.4 -0.2 0.2 

Xlr 



FIG. 3. The ratio of ReAo/Re^2 versus xlr for 771^(1 GeV) = 140 MeV. 



in the NDR scheme, where the calculations in the HV scheme are shown in parentheses 
Experimentally, the world average including NA31 |2T[, E731 KTeV |23| and NA48 
results is 



Re{e'/e) = (2.13 ± 0.46) x 10"^ 



(4i 



V. DISCUSSIONS 
A. Bag parameters Bi 

In Sec. Ill we have computed the non-perturbative parameters Bf'''^^ at /i = 1 GeV in 

(2) (2) 

NDR and HV schemes. Our results BI2 = 0.34 and Bg (q = 0.35 in the NDR scheme are 
smaller than the value 0.48 quoted in [|18[ for /i = 1.3 GeV. This is because we have taken 
into account isospin breaking contributions to the A/ = 3/2 amplitude so that IieA2 is 
enhanced by a factor of 1/(1 — Qm) [see Eq. ( [4.1| )]. Consequently, it is necessary to impose 
large nonfactorized effects and hence small i?|^to suppress A2. We observe from Tables II 
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0.08 0.1 0.12 0.14 0.16 0.18 0.2 

tris ( GeV ) 

FIG. 4. Direct CP violation e' je versus (in units of GeV) at the renormalization scale 
/i = 1 GeV, where the solid (dotted) curve is calculated in the NDR (HV) scheme and use of 
\a\{ytdyts) = 1-29 X 10"^ and Xlr = —0.1 has been made. The solid thick lines are the world 
average value for e'/e with one sigma errors. 



and III that the parameter Bg""* has the largest deviation from unity; it is equal to 8.7 in the 
NDR scheme. This is mainly because the ratio {Q2)o/ {Q2)o^^ = 1 — GXll [cf. Eq. (|3.7|)] 
is greatly enhanced by the nonfactorized effect, recalling that xll = —0.73 . Note that our 
results for b[^2 close to that obtained in the chiral quark model 0. 

There exist some nonperturbative calculations for B^^ and B^^l. Among them, lattice 
calculations are carried out at /x = 2 GeV and in the NDR scheme. However, the lattice 
results for B^fl are much more reliable and solid than B^\ Most approaches find B^'^ below 
unity and B^^ < Bq^\ while we obtain B^^ = 1.6 and B^^ ~ Bq'^ Since the radiative 
correction (as/4vr)[(f)87((57)2 + {'f)s8{Qs)2]/ {Q8)2'^^ is positive in our case [see Eq. (|3l^)], it 
is not possible to push B^ down below unity. 

As for the scheme dependence of Bi, it was arg ued in [EB] that sf^(HV) ^ 1.2 (NDR) 
and 5j^^(HV) ^ 1.2 sj^^ (NDR), whereas we find a very weak scheme dependence for Bq 
and 5f ^ but strong dependence for and 5}°'^^ ^^^(HV) = 1.7 (NDR), 5}°^(HV) = 
1.7S}°^(NDR) and 5f^(HV) = 1.8 (NDR). We have also studied the dependence of 
i?-parameters and found that only Bf'l 4 exhibit a significant dependence, while the other 
S-parameters are nearly independent. For example, we obtain B^'^ = 3.2 if ms(l GeV) = 
100 MeV and b[^^ = 2.5 if m,(l GeV) = 140 MeV (see Table II). Note that some other models 
predict a different rus behavior for i?-parameters. For example, Bq^^ is proportional to 
in the chiral quark model 0. 



B. X ^ vrvr amplitudes 

From Fig. 2 or Eq. (|4.6| ) we see that about (60-70)% of ReAg amplitude is accounted 
for in the present approach if ms(lGeV) lies in the range (125-175) MeV. Specifically, Qi, 
Q2 and penguin operators explain 66%, 18% and 14%, respectively, of the Aq amplitude for 
ms{l GeV) = 150 MeV. Hence, tree-level current x current operators account for around 85% 
of ReAg. However, contrary to [Q, we find that penguin-like diagrams induced by Qi, i.e.. 
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the penguin operators in Eq. ( |3.2| ), contributes only about 15% to RsAq. As conjectured 
in [0, the ly-exchange mechanism could provide an additional important enhancement of 
the Ao amplitude. Since the VT-exchange amplitude in charmed meson decay is comparable 



to the internal ly-emission one it is conceivable that in kaon physics the long-distance 
contribution to l^-exchange is as important as the external VT-emission amplitude. 

It is instructive to see how the predictions of ReAo and Rey42 amplitudes and the AJ = 
1/2 rule progress at various stages. In the absence of QCD corrections, we have 02 = |ai 
and 03 = 04 ■ ■ ■ = aio = under the vacuum insertion approximation. It follows from Eq. 
O that i 

° = = 0.9 (in absence of QCD corrections). (5.1) 
ReA2 4v2 

With the inclusion of lowest-order short-distance QCD corrections to the Wilson coefficients 
zi and Z2 evaluated at fi = 1 GeV, Aq/A2 is enhanced from the value of 0.9 to 2.0, and 
it becomes 2.3 if 771^(1 GeV) = 150 MeV and QCD penguin as well as electroweak penguin 
effects are included. This ratio is suppressed to 1.7 with the inclusion of the isospin-breaking 
effect, but it is increased again to the value of 2.0 in the presence of final-state interactions 
with 60 = 34.2° and 62 = -6.9°. At this point, we have Re^o = 7.7 x 10"^ GeV and 
IieA2 = 3.8 X 10~^ GeV. Comparing with the experimental values 

Re Ao = 3.323 x 10"^ GeV, Re A2 = 1.497 x 10"^ GeV, (5.2) 

we see that the conventional calculation based on the effective Hamiltonian and naive factor- 
ization predicts a too small A J = 1/2 amplitude by a factor of 4.3 and a too large A/ = 3/2 
amplitude by a factor of 2.5 . In short, it is a long way to go to achieve the A/ = 1/2 rule 
within the conventional approach. 

Replacing c]''^{fi) by the effective Wilson coefficients d^^, or equivalently replacing the LO 
Wilson coefficients by the NLO ones and including vertex-like and penguin-like corrections 
to four-quark operators, we obtain KbAq = 1.37 x 10~^GeV and R.eA2 = 3.3 x 10~®GeV 
and Rey4o/Re742 = 4.2. Finally, the inclusion of nonfactorized effects on hadronic matrix 
elements will enhance ReAo/ReA2 to the value of 13.8 with Re^o = 2.07 x 10~* GeV and 
Rey42 = 1.50 X 10~^GeV. To summarize, the enhancement of the ratio Rey4o/Rey42 is due to 
the cumulative effects of the short- distance Wilson coefficients, penguin operators, final-state 
interactions, nonfactorized effects due to soft-gluon exchange, and radiative corrections to 
the matrix elements of four-quark operators. Among them, the last two effects, which are 
usually not addressed in previous studies (in particular, the last one), play an essential role 
for explaining the bulk of the A/ = 1/2 rule. 



C. Direct CP violation e'/e 

From Fig. 4 or Eq. ([4.7|) we observe that, contrary to the case of A0/A2, the prediction of 
e'/e shows some scale dependence (see Fig. 4); roughly speaking, (£y£)NDR ~ 1-5 {e'/e)-^^. 
To understand this, we note that e'/e is dominated by be and bg terms (or y^^ and yf^) or 
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equivalently by the hadronic parameters Bq and (see Eqs. (^4.2| ), ( [4 .51) and Tables I-III). 
Moreover, direct CP violation involves a large cancellation between the dominant yf^ and 
yl^ terms. The scale dependence of the predicted e'/e is traced back to the scale dependence 
of the effective Wilson coefficient yf^ (see Table I). As mentioned before, formally yf^ should 
be scale independent to the order a^. It is thus not clear to us why ?/g^(NDR) and yf^(B.V) 
are not the same to the accuracy under consideration. Furthermore, the scale dependence of 
j/g^ is amplified by the strong cancellation between QCD penguin and electroweak penguin 
contributions, which makes it difficult to predict e'/e accurately. It appears to us that the 
different results of e'/ e in NDR and HV schemes can be regarded as the range of theoretical 
uncertainties. It is easily seen that a suppression of B^g'^ or an enhancement of Bq^^ will 
render e'/e larger; that is, a ratio of Bq'^ / B^^ greater than unity will help get a large e'/e. 
However, in our approach B^^ ~ Bq^^ = 1.5 and they are nearly scheme independent. We 
have also studied the dependence of e'/e on the nonfactorized effect xlr and found that it 
increases slowly with xlr (see Fig. 5), opposite to the case of Aq/A2. 

2 

1.75 

1 . 5 

„ 1.25 
O 

X 1 
0.7 5 
. 5 
0.25 

-0.6 -0.1 -0.2 0.2 0.4 0.6 0.8 

XLR 

FIG. 5. Direct CP violation e'/e versus xlr for lm{VtdVtl) = 1.29 x lO^"^ and ^^(1 GeV) = 140 
MeV, where the solid (dotted) curve is calculated in the NDR (HV) scheme. 

Since the AI = 1/2 rule and e'/e are both under-estimated theoretically, it is natural to 
ask if there exists a strong correlation between them. The two principal mechanisms respon- 
sible for the enhancement of Aq/A2 are the vertex-type as well as penguin-type corrections to 
the matrix elements of four-quark operators, and the nonfactorized effect due to soft-gluon 
exchange. Turning off these two effects by setting xll = Xlr = and y^^ — *• yf^, we find 
that e'/e does not get changed in a significant way. On the other hand, if a small strange 
quark mass is responsible for the remaining enhancement necessary for accommodating the 
data of A0/A2, it turns out that m,(l GeV) = 85 MeV and e'/e = (2.3 - 3.5) x 10"^ How- 
ever, this rris is too small even compared to the recent lattice result |^ which favors a lower 
strange quark mass: 771^(2 GeV) = (84 ± 7) MeV. 

VI. CONCLUSIONS 

The AI =1/2 rule and direct CP violation e'/e in kaon decays are studied within the 
framework of the effective Hamiltonian approach in conjunction with generalized factoriza- 
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tion for hadronic matrix elements. Our results are as follows. 

1. We identify two principal sources responsible for the enhancement of R,eAo/R,eA2: 
the vertex-type as well as penguin-type corrections to the matrix elements of four- 
quark operators, which render the physical amplitude renormalization scale and scheme 
independent, and nonfactorized effect due to soft-gluon exchange, which is needed to 
suppress the A/ = 3/2 X — > tttt amplitude. This approach is not only much simpler 
and logical than chiral loop calculations but also applicable to heavy meson decays. 

2. We obtain renormalization-scheme independent predictions for X — > tttt amplitudes 
and find ReAo/ReA2 = 13-15 if m,(lGeV) lies in the range (125-175) MeV. The 
tree-level current x current operators account for around 85% of RbAq. We conjecture 
that the VT-exchange mechanism may provide an additional important enhancement 
of the A/ — 1/2 amplitude. 

3. The bag parameters Bi, which are often employed to parametrize the scale and scheme 
dependence of hadronic matrix elements, are calculated in two different renormaliza- 
tion schemes by considering the vertex-like and penguin-like corrections to four-quark 
operators. It is found that B^^ ~ both of order 1.5 at = 1 GeV, are nearly 75 
scheme independent, whereas 53*^57 and Bj^^ show a sizable scheme dependence. Our 
results b[^1 = 0.34 and B^^Iq = 0.35 in the NDR scheme are smaller than the value 
quoted in the literature since we have taken into account isospin breaking contribu- 
tions to the A7 — 3/2 amplitude. As for the dependence of S-parameters on m^, only 
Si^3 4 exhibit a significant rus dependence, while the rest S-parameters are almost nis 
independent. 

4. Nonfactorizable contributions to the hadronic matrix elements of {V — A){y — A) four- 
quark operators are extracted from the measured K'^ — > tt+tt^ decay to be xll = 
—0.73, while the nonfactorized term for {V — A){V + A) operators lies in the range 
—0.45 < xlr < 0. We found that ReAo/ReA2 (s'/e) decreases (increases) slowly with 
Xlr- 

5. For direct CP violation, the prediction of e'/e is renormahzation scheme dependent 
owing to the scale dependence with the effective Wilson coefficient for reasons 
not clear to us. We obtain e'/e = (0.7 - 1.1) x 10"^ if m,(lGeV) = 150 MeV and 
e'/e — (1.0 — 1.6) X 10"^ if nis is as small as indicated by recent lattice results. 
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